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TECHNICAL NOTE 367 b 


THEORETICAL PRESSURE DISTRIBUTIONS FOR SOME SLENDER 
, WING-BODY COMBINATIONS AT ZERO LIFT 

By Paul F. Byrd 


SUMMARY 


Pressure distributions are calculated for some symmetrical wing-body 
combinations at zero lift. The theory of the calculations is based on 
the assumption of extremely slender wings and bodies and yields results 
for both subsonic and supersonic speeds. The examples considered are 
swept wings of constant chord mounted on bodies of nearly cylindrical 
form. 


Of particular interest is the effect of indenting the body on the 
distribution of pressure over the wing. When the indentation is such as 
to maintain a constant total area of the cross sections normal to the 
stream, the theoretical pressure disturbances remain small throughout the 
transonic range. With such indentation the isobars tend to remain smooth 
and nearly parallel to the sweep, of the wing surface. 


INTRODUCTION 


In several papers, important extensions to the Munk-Jones slender- 
body theory (refs. 1 and 2 ) for lifting wings and bodies have been made 
to include the theoretical effects of thickness on the aerodynamics of 
wings and wing-body combinations. Ward (ref. 3)> solving the linearized 
differential equation for the perturbation velocity potential by opera- 
tional methods, and employing asymptotic expansion of the solution, 
investigated the flow around bodies of general cross section at supersonic 
flight speeds. By a different procedure, similar results for a wing, 
body, or wing-body combination at subsonic speeds have been developed 
by Heaslet and Lomax (refs . 4 and 5 ) • An analysis for subsonic flow was 
also carried out independently by Adams and Sears (ref. 6 ) who, in addi- 
tion, made an extension for not-so -slender wings. Confining themselves 
to wings at zero angle of attack, Keune (ref. 7) and Oswatitsch and Keune 
(ref. 8 ) have recently obtained a slender-body theory that is slightly 
different from those of references 4, 5, and 6 . In a more recent report. 


Supersedes NACA RM A5^J07 by Paul F. Byrd, 195^- 
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Harder and Klunker (ref. 9 ) have applied the basic ideas of the slender- 
body approximation to the nonlinear transonic equation for the velocity 
potential. 

The principal object of the present investigation is to apply the 
general method of reference 4 in calculating the pressure distribution 
for some special cases of nonlifting slender wing-body combinations in 
subsonic and supersonic flow. The wing of the combinations is swept back 
and has a symmetrical section with rounded leading edges. Determination 
is made of the pressure for the wing alone and for cases when the wing is 
mounted on a circular cylinder or combined with a body indented such that 
the axial variation of cross-sectional area of the combination is constant. 
The effects of Mach number and sweep angle are included in the results 
presented. 


LIST OF IMPORTANT SYMBOLS 


^o 

c o 

C P 


(V 

Cn 


L o 

m 

M 

r 

Re 

Ro 

R(x) 

sgn(x-g) 

so 

s(x) 

S(x) 


value of x at which s Q = s(x) 
root chord 


pressure coefficient. 


-2u 
U n 


(Vy 2 + 


pressure coefficient on indented wing-body combination 
over -all length of the wing 
slope of wing leading edge (See sketch (f).) 
free -stream Mach number 
polar distance in y, z plane 
real part of a complex quantity 
radius of cylindrical portion of body 
radius of indented body of revolution 
sign of (x-§) 
maximum value of s(x) 
local semispan 

local cross-sectional area of wing alone 



MCA TN 3674 


3 


t(x) 

u 


Uo 


V 


v r 

w 


z 


P 

0 

A(x,y) 

<P 

£ 


B 

I 


sub 

sup 

u 


W 


spanwise distance from x axis to wing trailing edge 
perturbation velocity in the x direction 
velocity of free stream 

perturbation velocity in the y direction 

radial component of velocity in yz plane 

perturbation velocity in the z direction 

Cartesian coordinates (x downstream, y to starboard, z 
upward) 


- 1 | 


polar angle in yz plane 

slope of wing surface in x direction 

perturbation velocity potential 

part of potential satisfying cpyy + cp zz = 0 

complex variable (y + Iz) 

complex variable (y* + izi) 

constant related to t q 

maximum thickness of wing section 


Subscripts 


body 

lower surface of the wing ( z = 0 plane) 

subsonic 

supersonic 

upper surface of the wing (z = 0 plane) 
wing 
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ANALYSIS 


The thickness distribution of a symmetrical- wing is prescribed. 

Under the restriction that the thickness is small and that the configura- 
tion is slender, formulas to be applied later (see "Applications" section) 
will now be briefly presented for determining the pressure coefficients 
for particular nonlifting wing-body combinations. Equations are first 
given for the wing without a body, and are then modified for cases when 
the wing is mounted on a circular cylinder or on an indented body of 
revolution. 


Nonlifting Wing Alone 


The differential equation and boundary conditions .- Expressed in 
terms of the perturbation potential cp(x,y, z), the basic linearized par- 
tial differential equation for subsonic as well as supersonic flow is the 
familiar Prandtl-Glauert equation 


(! - Mq 2 )^ + CPyy + CP ZZ = 0 


( 1 ) 


where M 0 is the free -stream Mach number. If the surface of a wing 
z(x,y) is given, solutions to the differential equation must satisfy the 
boundary condition that the flow is parallel to the wing surface. When 
the wing is thin, it is sufficient to satisfy this requirement in the 
plane z = 0. Analytically, the expression of the condition is 



v u(x,y) 



UoM x >y) 


( 2 ) 


where Uo is the free-stream velocity and w u (x,y) is the vertical 
induced velocity on the upper side of the z = 0 plane. 

Velocity potential .- When the flow is supersonic, the formula for 
the perturbation potential subject to the boundary condition (2) is known 
to be (e.g., see ref. 10) 


<p(x,y,z) = - yy\i(s>ii) 

T 


arccosh 


* - i 


P-s/ (y - ti) : 


dT] (3) 


+ z* 
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where the region of integration T is the portion of the plan form lying 
within the Mach forecone from the point x,y, z. In the case of subsonic 
flow fields, the solution of the differential equation (l) may he expressed 
in the form 


<P(x,y,z) = - ST ^ J'J' arcsinh 

T 


x - 6 


Ji: 


- Tl ) 2 + z 2 


<11 dr] (4) 


with T now extending over the entire plan form. These two solutions 
yield the potential due to a distribution of sources of strength propor- 
tional to the slope of the wing surface \i(x,y) . 


When the wing is slender, that is, if 


3 J(y - t) 2 + z 2 

x - i 


is considered 


very small, further approximations to the linearized potential ( 3 ) and ( 4 ) 
can he readily made. In this event, one may employ the approximate 
relations 


arccosh 


x - i 


3 Ji y - T]) 2 + 


£ In 


2(x - |) 

3 J (y - ti) 2 + 


> (5) 


arcsinh 


x - 5 


3«/(y - Tj) 2 + 


- af sgn(x - £)ln 


g|* - 6| 

pVty - i ) 2 + z 2 J 


where the symbol sgn(x - §) means that the sign of (x - |) is to he 
taken. 


Consider now a thin pointed wing 
of symmetrical section with straight 
or sweptforward trailing edges as is 
shown in sketch (a) . Use of the 
relations (5) in equations ( 3 ) and (4) 
then gives the result (ref. 4) 


9 (x,y,z) =<p 2 (x,y,z) + g(x) 


(6) 



Sketch (a) 
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where 


Uo ^?3( x ) 

<P 2 (x,y,z) = / A u (x,T])ln[(y - T }) 2 + z 2 ]dq 

s 2 (x) 


and 


g(x) 


Uo 

bx 



2(x - 0 
P 


d| 


(7) 


(8a) 


for supersonic flow, and 

Uo b r° 3 lx - 1 1 

g(x) = - ^ sgn(x - |)S*(|)ln - dg (8b) 

o 

for subsonic flow. The function S*(x), the derivative of the cross- 
sectional area of the wing in a. yz plane, is found from 


S*(x) = 2 


_d 

dx 


s s (x) 

/ Z u (x,.Tl)dTl 


(9) 


s„(x) 


It is seen from equation (6) that the slender-body approximation to 
the linearized supersonic and subsonic potentials consists of two parts. 
The first part (<P 2 ) in each case is independent of Mach number and is a 
harmonic function in the transverse plane; that is, it satisfies the two- 
dimensional Laplace equation 


'Pyy + tyzz - 0 


(10) 


The second part (g) depends on the cross-sectional area of the wing and 
is a function of x and f3 only. 

Inspection of equations (6) and (8) shows that the value of the 
potential for a particular wing at some Mach number Mi can be written 
in terms of the potential given at another Mach number M 0 . Thus 
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^PPi 9su Pp o + 2^ S ’ (x)ln ^ 


U 0 8t 

cp 8ubp i = cp sub Po + ST s '( x ) ln 


(11a) 


with Px - y|Mx 2 - l| and p Q = </" |^o 2 - l| . The supersonic and subsonic 
potentials may also be related. If one has already found an expression 
for the supersonic potential, he can then obtain the subsonic potential 
by means of the equation 


^sup + 2n S’(x)ln 


+ S/ 


S»(l) 
x - S 


(lib) 


Certain symmetrical wings whose trailing edges are swept back as in 
sketch (b) can also be treated by the simplified theory, provided the 
chordwise variation in shape of the 
cross sections is sufficiently smooth \Uo 
and gradual that the assumption of 
two-dimensional transverse flow may 

reasonably be applied. Since the T yTC — 

wing is at zero angle of attack and / \ 

the flow is symmetrical about the fy~~ s M/ \s~y=s(x) 

xy plane, the source distribution C ° \ 

A u (x,y) in the plane of the wing in / / \ 

region W between the axis and the / ^ \ a ° 

trailing edge is set equal to zero. \ 

If the limits of integration are / S' \ 

properly adjusted, equation (6) will / 

then formally still apply and yield Ls _ _ 

expressions for the potential in the * -So ^ 

two regions 1 and 2 of the sketch. 

Thus for region 1, in the z = 0 plane, X 

one obtains Sketch (b) 


-y=s(x) 




u o r s 

'Psup _ ~ J Ajx^lnjy 2 - Tl 2 |drj - 


-y=t(x) 


^ J Sx l (|)ln | (x - £)d|, (x<c Q ) (12a) 
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and in region 2 


<P 


sup 


u o r s u o 8 r° 

= If J A u (x,Tj)ln |y 2 - n 2 !cLTi - ^ ^ Si^gJinCx - |)d| - 


H j r S2 . ( I ) ln| ( x- IHS, (x > o 0 ) 


(12b) 


where 


s(x) 


s i'(x) = 4 ^ J' z u (x j> T])dii, 0 < x < c Q 


s(x) 


S 2 *(x) =4“^ ZuUjTlJdl), C Q < X < J 
t(x) 


(13) 


The subsonic solution is obtained from the supersonic solution by using 
equation (lib). 

Pressure coefficient .- After determining the potential for a wing 
from the equations in the foregoing section, the pressure on the surface 
of the wing is found by differentiation. The pressure coefficient is 
related to the perturbation velocities by the equation (ref. 4) 


c p (x,y,0) 


_ 2 _ 85 1 /dqA 2 1 

_Uo 8x + u7 V W + U 0 2 



(l4a) 


which is an approximation to the complete Bernoulli equation consistent 
with both the linearized differential equation (l) and the assumption of 
slenderness. For planar problems, further simplification achieved by 
neglecting the nonlinear terms (dcp/dy) 2 and (8<p/dz) 2 yields satisfactory 
estimates for the coefficient. 
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Thus 


c p (x,y,0) - - (|f) < lto > 

z=0 

where <P is obtained from equations (llh) and (12). 

It should he pointed out, however, that equations (lib), (12), and 
(14) will not furnish realistic results for the pressures 1 for a wing in 
either subsonic or supersonic flow unless certain restrictions are imposed 
on the gradient of cross-sectional area and its derivative S"(x) . Con- 
sider, for example, the pressure found along the line AB where regions 
1 and 2 in sketch (b) join, 
tions (12) and (l4b) that at 
sure coefficient given by 


Ac p = c Pi - cp 2 = ~ 


4 


In supersonic flow, one finds from equa- 
this line there is a jump Ac p in the pres- 


lim 
e— > 0 


|aS"( c 0 


)ln 


2e 


P's/ s2(c 0 ) - y 2 


dt 

dx 


A u (x,t) 


In — 
C n +6 ^ 


( 15 ) 


where 


AS"(c 0 ) = Si"(c 0 ) - S 2 "( c G ) 


and where the usual assumption that S’(0) = S*(c 0 ) = 0 has been made. 
Evidently, along the line AB a logarithmic singularity will occur in 
the pressure distribution in going from region 1 to region 2 if 
Si"(co) / S2 ** ( Co) • (The infinity would of course be higher than loga- 
rithmic if AS"(c 0 ) is singular.) The formulas for the pressure coef- 
ficients presented in this report and in reference 4 are therefore good 
only ^or cases where the plan form and the slope of the wing surface are 
sufficiently smooth so that there are no abrupt changes in either S’ or 
S". (This holds true also for the formulas given in references 3> 6, and 
7. The restriction on S", however, may be somewhat relaxed in employing 
the slender -body theory of reference 8.) 

1 Although the theory may give spurious infinite pressures on certain 
portions of the wing, the results obtained for the wave drag by integrat- 
ing the product of pressure and surface slope over the wing may be finite 
and reasonable . 
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Even if there is no jump in the pressure coefficient at the line AB, 
the slope of the pressure curve may not he continuous. It can be shown, 
for instance, that a singularity will occur in the slope of the curve for 
the pressure coefficient in case S n '(c 0 ) is discontinuous. 


Wing on a Cylindrical Body 

The equations presented in the foregoing sections for the wing alone 
will now be modified to yield formulas for calculating the pressure dis- 
tribution on a combination composed of a 
| Uo symmetrical wing z u (x,y) = -z z (x,y) 

mounted on an infinite circular cylinder 
having radius R 0 (see sketch (c)). The 
surface of the wing chosen here will also 
be considered symmetrical about the xz 
plane, but the method applies equally well 
if this is not the case. The procedure 
followed is essentially the same if the 
fuselage is any body of revolution instead 
of a circular cylinder. 

In studying such combinations, it is 
usually convenient to introduce a second 
coordinate system. Let the yz plane be 
represented by a complex variable 

£ = y + iz = re 1 *-® (l6) 

and then consider a £1 plane 
£1 = yi + izi = rxe 10 ! 

obtained from the £ plane by the 
Joukowski transformation 
Sketch (c) 2 

q - i * \ (17) 

The transformation (17) maps the £ plane onto the Ci plane so 
that a circle representing a section of the body in the £ plane is 
mapped onto a portion of the real axis in the £3. plane, while the part 
of the real axis outside the circle is transformed into an adjoining part 
of the real axis of the Ci plane. (See sketch (c).) It can easily be 
shown from the equation that the geometric relations 
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Ri 

yi 

yi 

Sl 


2Rr 


y + 


Rr 


y ' 

2R 0 cos 9, 

R 0 2 

s + — ~— 




yi 


2 > R1 2 


yi 2 < R1 2 




(18) 


hold for zi equal to zero. 

The technique employed is to transform the boundary conditions from 
the £ plane to the £ 3 . Pl an e and to find a solution in the latter 
plane. The solution is then transformed bad; to the physical plane for 
the completion of the problem. 

Perturbation velocity potential at the wing surface . - A consequence 
of the conformal transformation is that the complex velocities in the 
two planes £ and are related by the equation 


<i£, 

v - iw = (vi - iwi) -£jr- 
or, in polar coordinates (ref. 11 ), 


(19) 


v = v x 


W = Wi 


1 - 1 5 ° 

-L ' r 


R, 


cos 29 


1 - [ cos 29 


1 A f R oY 

'J + Wl 

,1 w 

’J - V 1 -J 


sin 29 


sin 29 


> (19a) 


v r = 


Vi cos 9 + Wi sin 0 


1 - 


R N 2 
n o 


+ 2wi (^-^) sin 9 


From these equations it follows that the boundary conditions in the Zi = 
plane are 
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wi 


v r = U 0 (dEo/dx) 
2 sin 9 2 sin 9 


yi 


< R1 2 


> ( 20 ) 


Vi 


W 

l - Ro 2 /^ 


R1 2 < yi 2 < si 2 


with w related to the streamwise slope of the wing surface hy equa- 
tion (2) . 

The two -dimensional solution Cfe in the £1 plane at zx = 0 is, 
from "equation (7)> 


8 1 

qp 2 (x,yx, 0 ) wj.(x,T|x)ln|yi -T|il d ni 


( 21 ) 


Putting now 


“Si 


tii = n + 


*o 2 


and using relations (18) and the houndary conditions (20), one finally 
has for equation (21) 


tt s ( x ) 

U 0 p 

<P 2 (x,y) = IT A u (x,q)ln 


Rr 


(y 2 - ti 2 ) (yV - R 0 4 ) 
y 2 ^ 2 


dq 


(22) 


If one of the functions g given by formulas (8) is then added to this 
equation for cp 2 , the velocity potential on the wing may be obtained for 
either supersonic or subsonic flow. Thus, for the supersonic perturbation 
potential on the upper surface, there results 

u 

<psup = irf 

Ro 


(y 2 - ii 2 )(y 2 ii 2 - Ro 4 ) 
y 2 ^ 2 


dq 


H2 r s»(D 

2 « dxg 


in | (x - g)d£ 


(23) 


and for the subsonic potential 
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9s ub 


p) 

= ir j V(x,i 


Ti s )(y 2 ii 2 - R 0 4 ) 

- p - 2 dT ) 


L o 

[ sgn(x - |)S»(£)ln f |x - ||d£ 
4-it OX J B 


Pressure coefficient at the wing surface .- The pressure coefficients 
on the surface of the wing of the combination are obtained by differenti- 
ating equations ( 23 ) and (24) with respect to x and then employing 
relation 2 ( l4b) . Thus 


2 *r\, s (y 2 - n 2 ) (yV - V) 

11 2>xJ 3^2 


(x - £)a6 


(25a) 


^Psub 


2 a ° r (x) (y 2 - n s ) (*V - Ho 4 ) 

■«a -J Jk “ (x ’’ l)ln p? d,1 + 


oo 

■V p A 

Sx 2 J s e n ( x “ Us l (£)in | |x | d| 


(25b) 


or, with the aid of equation (lib). 


Psub 1 s r°s»(£) 


1 r. fl / \-, ^SUb 1 C> P 

= =p Bup - * < x)ln u ' 2* a x J 


As in the case of planar problems, the squared terms in the pres- 
sure relation can be neglected in considering a combination whose body 
is a circular cylinder. (See ref. 4.) 


Ik 
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The coefficient for a particular wing at some Mach number Mi, written 
in terms of the coefficient at another Mach number M 0 , is 


and 


C P 

SUp Pr 

Cp subp 1 


°p 

SU %o 


+ 1 s "^ | 

+ i S"(x)ln ^2 


(27a) 

(27b) 


where Pi = J~\ Mi 2 - l| and P D = J (Mq 2 - 1 1 . 


It should be noted that 
valid for all values of R n . 


\Uo 



formulas (25a) and (25b) are not uniformly 
If the radius R 0 is equal to zero, the 
■ equations, reduce immediately to those 
given previously for the wing without 
a body. As the radius approaches 
infinity, however, one finds that for- 
mulas (25a) and (25b) yield a value of 
the pressure coefficient that differ 
from the results for the wing alone by 

a term equal to — * — - 
the equations for values of 


it ln 2 


Sketch (d) 


Employing 
R 0 that 

are large violates the assumption of 
slenderness according to the approximate 
theory used herein. 

When the wing mounted on the com- 
bination is swept back, as shown in 
sketch (d), the slope A u (x,t|) is taken 
equal to zero in the gap between the 
trailing edge and the body. The for- 
mulas obtained in regions 1, 2, and 3 
for the supersonic pressure coefficient 
on the upper surface of the wing are 
then as follows: 


Region 1: 


2 d P ^ 
°Pl = -*5 xf 


(y 2 - ti 2 )(y 2 n 2 - R 0 4 ) . 

yV 


dr| + 


d 2 


1 

51 Sx 2 


J' Si 1 ( i ) ln j| (x - 5)dS, 


(0 < x < Co) 
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Region 2: 


s(x) 

i tf Au(x ’' l)ln 

t(x) 


(y 2 - •n 2 )(y 2 n 2 - Ro 4 ) a 
— 6 xi + 


« I?/ Sl ’ (l)ln(x ‘ ^ + « I?/ s 2 , (^ ln | U - ^ 


(c 0 <x<b 0 ) (28b) 


8 P° 

’Ps = - I ^ / *uU,Ti)ln 

t(x) 


(y 2 - Ti 2 )(y 2 T] 2 - R 0 4 ) „ 

dr, + 


Cq dq 

Si‘(^)m(x - |)dS + l ^ f S 2 «(i)ln(x - s)d£ 


x 

| /s a «(i)m | ( X - i)ai, (b 0 '<x<2 0 ) 


(28c) 


where b Q is the value of x at which s G = s(x) . Use of the relation 
(26) in conjunction with the above three formulas will furnish equations' 
for the pressure coefficient in subsonic flow in the various regions. 


Wing on an Indented Body 


In the previous formulas for the pressure distribution attributable 
to thickness, the coefficient becomes infinite 3 when the Mach number 

approaches unity because the formulas contain a term involving S T (x)ln — . 

P 

It is thus possible to construct slender configurations which will give 
a theoretically finite pressure coefficient even at the speed of sound. 
Under the assumptions of slender-body theory, this may be achieved if the 
gradient of cross-sectional area of the configuration in a yz plane 
vanishes identically. Combinations constructed in this manner on an 


a 

This of course is not a property of slender-body theory alone. 

Except in particular cases, steady-state linearized theory also in general 
ields infinite pressures at P = 0. 
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inf inltely long fuselage are slightly indented in the vicinity of the 
wing and possess, according to the theory, the important property of 
having zero wave drag. 


Velocity potential at the tying surface .- Let the surface 
z u (x,y) = -z i(x,y) of a symmetrical wing be specified and assume that 
the indented fuselage is a body of revolution which deviates Blightly 
from a basic circular cylinder of radius Rq . The cross-sectional area 
S(x) of the exposed wing is then given by 


\Uo 



or 


S(x) 


= k 


(x) 

z u (x,Tl)dT| 


R(x) 


(29) 


where R is the radius of the indented 
body. (See sketch (e) .) In order for 
the streamwise gradient of cross- 
sectional area of the entire combina- 
tion to be zero, the relation 


«R 2 ( x) = itR 0 2 - S(x) , S(x) < < rtR 0 2 

(30a) 


= -S'(x) 


(30b) ’ 


must hold. Since the quantity S(x) is known from formula (29), equa- 
tion (30a) can thus he solved for the radius R of the body. 

Now the perturbation potential in the plane at zi = 0 has the 

form 


or 


<P 2 (x,yi, 



s i 

■wi(x,Tii)ln | yi 


-Si 


- Til | dqi 


(31a) 


<P 2 (x,yi,0) = v 1 (x,T]i)ln|y 1 - qijdqi + 

-n r ! 


, s i 


^i ( x, t]i ) la |yi 2 - t)i 2 | dqi 

(31b) 
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with 


R 2 

si = s + — , 


R 2 

yi = y + -y- > ‘ r i = 2R 


(32) 


The term corresponding t6 the g(x) of equations (8), which is to be added 
to this <p 2 , is zero since there is no change in the cross-sectional area 
of the specified shape of the yz plane. 


From the last relation in equation (19a) it follows (upon putting R 
in place of Rq) that in the interval- tji 2 < ri 2 tne boundary condition 
may be expressed 


v r (dR/dx)U 0 

2 8in e = 2^/^^TW^IP ; 


(33) 


For very thin w i ng s , and for deviation Rq - R of the same order aB the 
wing thickness, the condition in the interval ri 2 < % 2 < B ! 2 is approxi- 
mately 


v _ Up(dz u/dx) 

1 - R 2 /^ 2 1 - R 2 /ti 2 


(3b) 


Making use of equation (33 ) , 
the right in equation (31b) as 


one can then write the first term on 


r x 

wjx^jln | 7l - Tijdr^ = UoRg lay, y 2 > R 2 ‘ (35) 

where, from equation (30b), the quantity R(dR/dx) may be replaced by 
-l/2jt times the gradient of cross sectional area S’ (x) of the wing . 

In the second integral in equation (31b), set 


Tli = n 


+ 


Rf 

3 


and use relations (3^). The final form for the velocity potential at the 
upper surface of the wing thus becomes 

U o 

<P(x,y,0) = " ^ S’ (x) In y + 


u o 

— / A u(x,T])ln 

R 


(y 2 - r 1 2 )(y 2 T 1 2 - R 4 ) 


y 2 T) 2 


dii, 


y > R 


(36) 
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Formula ( 36 ) for the velocity potential for a slender combination 
whose body is indented according to the area rule (ref. 12) is evidently 
independent of Mach number and therefore holds (within the assumptions 
of the simplified theory) for subsonic, transonic, and supersonic speeds. 
If the limits of integration are adjusted, the equation may also be used 
for combinations with sveptback wings . 

Pressure coefficient at the wing surface .- On the surface of the wing 
the pressure coefficient can be found from relation (l4a), that is. 


2 dcp 1 /clcpN 2 ‘ 

%jr |_U 0 chc + U Q 2 \~dyj J 


where the potential <p is obtained from equation ( 36 ) . 


APPLICATIONS 


The formulas which were given the preceeding part of this report will 
now be applied for the purpose of performing detailed calculations of the 
pressure distribution for some particular nonlifting combinations having 
a symmetrical wing with constant chord. 


Wing Mounted on a Body of Circular Cross Section 


Pressure coefficient on the surface of the wing .- The wing of the 
combination 4 considered here is mounted on a circular cylinder and has 
an upper surface defined by 

z u (x,y) = T(c 0 m - s + y)v/(s - y)(c Q m - s + y) (3^ 


s = mx + R, 


zu( x >y) = z u (x,-y) = -z j(x,y) = -z z (x,-y) 


8t d 

T = 

3^0 a/ 3 - 


■^This combination will be referred to as the basic combination. 
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the constant r Q being the maximum thickness . The profile resembles a 
Joukowski section and is the same for all values of y. On the st ar board 
side of the wing , the slope of the surface is given by the relation 


5z u 

8x 


A u (x,y) 


w u (x,y) 

U 0 



- 4s + 4y) 


/ 


c Q m - s + y 

s - y 


(40) 


From this equation it is seen that the slope possesses a square-root 
infinity at the leading edge and is zero at the trailing edge. 

The gradient of cross-sectional area of the wing is 
Si’(x) = 4Tm 3 (c Q - x)*»/x( c Q - x) = 4m z^UjRq), (0 < x < c 0 ) (4la) 

in region 1 in sketch (f) and is- 

S2 ’(x) =0, c D < x < (s Q - Ro)/m (4lb) 

in region 2; for region J,, 

^3 (^) = “4m Zjj(x^Sq)^ (sq ~ Rq) /m ^ x ^ (bq - Rq + CQm) /m 


Since S" is not continuous across 
the line AB, the approximate theory 
employed here will give unrealistic 
results in the regions of the tips. 

In calculating the pressures for the 
wing-body combination, attention will 
therefore be confined to regions 1 
and 2 only. The calculations to be 
given in subsonic flow will be based 
on the assumption that the tips are 
located far downstream and have no 
effect on the other two regions . 

Use of the equations (40) and 
(4l) in formulas (28a) and (28b) yields 
the following final results for the 
supersonic pressure coefficient on the 
upper surface of the wing for regions 
1 and 2: 


I 


U 0 


(4lb) 



Sketch (f) 
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Region 1: 


'Psup 


= 2— l6R n tan -1 / — - ■ — ' + 2m(co - 4x) j 
it L 0 v c D - x v 


c n - x 8xy(c 0 - x) 

-2= In — — + 


PcoCy 2 - R 0 2 ) 


F(x,y) + F(x,-y) - 2F(x,0) + F ( x ^r) + f(x, - 

(0 < x < c Q ) (42a) 


Region 2: 


°V 


sup 


= n^T 


8r, 


/ x - c n ,, „ /s - c 0 m 

o + 2m(4x - c Q ) J — + 2(4s - c 0 m)y- 


s 


/ s - c 0 m + y / 

(c Q m - 4s - %-y + (^ c o^ 4s + 4 


, RcA / s - c Q m + Ro2/y 

c_m - 4s - 4 ~rr- ) / ~ 

0 J JJ B + R^/y 


Bp \ /s ~ c Q m - Rp g /y 
y s - R 0 g /y 


(c Q < x, s - Com < y < x; s = mx + R 0 ) 

(42b) 


where 

y s - c 0 ia + 

r~ 

s + ij 


[§- 


tan 


_i (c Q - x)(2s + n) - (tjx + RqCq) 
2 *Vx(c 0 - x)(s + T])(s - c Q m + Tl) 


] 


(43a) 


when c m - s - ij < 0, and 


F(x,tl) = (c Q m - 4s - 4H)/ ° s ‘ “ 


- B - T] 
*1 


In 


c O^ T l + R o) 

(c Q - x) (2s + Tl) - (tlx + RqCq) + Vx(c 0 - x)(s + ti)(c 0 tn - s - tl) 

(43b) 


when Cota - s - t| > 0 
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It is seen from equation (2 6 ) that the formulas for the pressure 
coefficient in subsonic flow can be obtained by adding two terms to the 
above equations for the supersonic pressure coefficient. Since the value 
of S’ in region 2 is identically zero, the terms to be added for region 
1 become 

c 0 

" hi 'Sxf x - ~ g ^ = TmS ^ x - 3c 0 ], (x < Co) (44a) 

o 


and for region 2 


~L 

2a 



rm 


4x - 3c 0 + (c Q 



(x > Co) 

(44b) 


The subsonic pressure coefficient on the upper surface of the wing can 
therefore be written in region 1 as 


-sub 


m^ 

°Psup + it 


it(4x - 3c Q ) + 2(4x - 



(45a) 


and in region 2 as 


Cp = c_ + m 3 T| 


sub P 


sup 


4x 


- 3c 0 + (c Q - 4 x)J - ~ C °- 


(45b) 


where c^ 

-^sup 


is given by equations (42) . 


Formulas (42) and (45), which give the pressure coefficients on the 
wing of the combination indicated in sketch (f), will be plotted and the 
results discussed in a later section. An obvious result noted now is 
that in region 2 the formulas do not depend on Mach number for either 
supersonic or subsonic flow, but that the subsonic coefficient in this 
region is always greater than the supersonic coefficient because the 
inequality 


m 3 T 


4x - 3c 0 


+ (c Q 



(46) 


holds for all values x ^ c 0 . 
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As a check, let us now consider the asymptotic behavior of equa- 
tions ( 38 ), (42b), and (45b) far outboard along the wing. For this pur 
pose, we introduce a change in coordinate system by means of the trans- 
formation 


x = x n sin 0 + ynsin 0 
y - Ro = -x n eos © + y n sin 0 


U 0 


and 



(See sketch (g).) Setting tan 0 = 
one can then write 


x = 


mx n + y n 


y - Ro =■ 


\l 1 + m 2 
Wu - x n 


*n 


-Jl + m? 
mx - (y -Ro) 




\J 1 + m £ 


= 


x + m(y - R 0 ) 


-Jl + m 2 


J 




= c 0 m/Jl + m 2 


% = m/Vl + m? 

Uq = U 0 m/s/l + m 2 


p n = l/^/l + m 2 


(i 


Equation ( 38 ) thus becomes 


8 T o , 3/s, n 1/2 

z - = — F^-~ ( c n “ Xn) (xn) 


n 
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which agrees with the approximate equation employed in two-dimensional 
section theory for a thin Joukowski base profile of thickness ratio 

— — . (See ref. 13 .) The pressure coefficient for such a two- 

3*/3c n 

dimensional airfoil is 


2 dtp 

^ " Ua Sx n " 

or, witten in terms of the xy 


8 tq 

3 'V3cn Pn 




system. 


°P “ " 


8 t 0 

3 4 3c 0 2 


[3c 0 m - 4(mx - y + K 0 ) ] 


(49a) 


(49b) 


This equation is in agreement with the asymptotic expression of the pres- 
sure coefficient for region 2 obtained far outboard along the wing from 
either equation (42b) or (45b) . 

Pressure coefficient on the body .- The formula for determining the 
pressure coefficient on the surface of the body is (ref. 4) 


% 


2 Sep 1 
|_U 0 cix U 0 2R £ 



( 50 a) 


When the body is a circular cylinder, this relation reduces to 



( 50 b) 


It is apparent that application of equation (50b) requires a knowledge 
of the value of the potential 9 in space, so that formulas ( 28 ), which 
hold for z = 0, must be "analytically continued." One uay that will in 
effect accomplish this is as follows: 

Set 


y a ■ d i® 

i = y + iz = R 0 e 


and form the function 

®(x,£, 0 ) = 9( x ,£,0) + it(x,£, 0 ) 


(51) 



NACA TN 367^ 


where cp(x,{;,0) is taken from equations (23) and (24), with the function 
i|r defined by 5 


s(x) 

t(x,^,0) = J w u (x,q)dTi 


The constructed function (51) will obviously satisfy Laplace's equation 
in the yz plane; and it can be verified without difficulty that the 
boundary condition v r = 0 is also satisfied. The coefficient on the 
body can then be calculated from the formula 


2 w f ^ 
°p " u« Re V ^ 


where the symbol Re means that the real part is to be taken. 

For region 1, the final result obtained for the pressure coefficient 
on the body in supersonic flow for the basic combination considered in 
this part of the report is 

^sup ■ UT + F(*,R 0 e") + 


F(x,R o e' 10 ) + F(x,-Roe 10 ) + F(x,-R o e' 10 ) + 


, , % / c D - x , 8x(c 0 - x)e 

Re 2m(c 0 - kx) J — ^ —TZ — 

pc 0 (e^ - 1)R 0 


iitCcQm - 4s + 4R Q e i ) 


'c D m - s + i^e 
s - Roe 10 


U < c 0 ) 
(54a) 


where the function F(x,q) is given by equations (43) . The result found 
for the supersonic pressure coefficient in region 2 can be put in the fprm 


s An arbitrary real function E(x) may, in general, be added to the 
right side of equation (52). Such a function, however, would in no way 
affect the real part of B$/dx, since iE is purely ima 


E tsas a 
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°P 


= m?T 


sup 


8R 0 + 2m(4x - c G ) J - ^ C ° - + 2(4s - c Q m) I — + 


>/ 


±e 


(c 0 m - 4s - 4Roe °) 


i0 b + Rpe - c 0 m 


s + R 0 e 


10 


(com - 4s - 4R 0 e i0 ) 


I s + Rn^ ^ " C n lll 


s + Roe 


-10 


(c Q m - 4s + 4Roe i0 ) l~— ° ~ + 


s - R 0 e 


10 


(c Q m - 4s + 4R 0 e -:i -®) 


Ib - R o e“ i0 - c Q m 
s - R 0 e“^® 


(x > c 0 ) (5^) 


Corresponding formulas for the pressure coefficient on the body in sub- 
sonic flow are obtained from the above two equations by adding to them 
the terms on the right of equations (44a) and (44b), respectively. For 
any fixed value of y in region 2, relations (45b) and (46) furnish the 
inequality 


'T’sub > ^Psup 


(55) 


which, as mentioned before, is also true for the coefficients on the wing. 


Wing Mounted on an Indented Body 


The constant-chord wing whose thickness distribution on the starboard 
side is given by 

z u (x,y) = T(e 0 m - s + yW (s - y)(c Q m - s + y) , R(x) ^ y < s j 


s = mx + R 0 

t 8t q 

31a 2 Cq 2 s[3 


>(56) 
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will now be combined with an indented body of revolution which deviates 
slightly from a basic cylinder in the manner that makes the local cross- 

sectional area of the configuration 


\U a 



a constant. The wing has the same 
profile along EE as the one in the 
previous example. (See sketch (h).) 
Since the gradient of cross-sectional 
area of the combination is set equal 
to zero , the wave drag, according to 
the slender-body theory, is also zero. 

The cross-sectional area of the 
w in g in region 2 is equal to the con- 
stant 


s 

Sg(x) = 4 J - 


(57) 


s-c 0 m 


end the radius of the indented body 
in this region is therefore the con- 
stant 


“a 


- “ a o J 

V *R 2 V 


Cq 3 ™ 3 

too 2 


(58) 


Corresponding equations for region 1 are 


s • f 

Si(x) = h J P z u (x,T 1 )dTi = | 3 m s c 0 


c -2 c 0 m - 2s + 2Ri , 
a cos i + 


Co m 


2[2(s - Rx) (7 c 0 m - 4s + 4 Rx) - 3c 0 2 m 2 ] N /7T^^Kcom^TTR7)j-(59) 


and 


Rx(x) 


= Ro 



0 — Ri ^ Ro 


(60) 
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Equation ( 60 ) , however, does not give an explicit value of the radius 
in region 1 because the function S x itself involves R x . Moreover, 
the function Si is a transcendental function of R^., so that equa- 
tion ( 60 ) can only be solved by graphical or other methods of approxima- 
tion. By means of the mean- value theorem, we can write the equation in 
the form 

Sjl(x) = S xi (x) + 4z u (x,d)(R Q - R x ) (6l) 

where d is a certain value in the interval R x < d < R Q , and 


Sn(x) 


m^T 

T2 


2(l4xc 0 - 8x^ - 3co z )4 x(c 0 - x) + 3c 0 s cos -1 



( 62 ) 


is the cross-sectional area of the wing in region 1 neglecting the addi- 
tional area exposed by the indentation. Since the deviation R 0 - R x is 
assumed to be of the same order as the wing thickness, neglect of the 
term 4z u (x,d)(R Q - R x ) will evidently introduce only an error of the 
second order in thickness. For very thin wings, the radius E x may thus 
be approximated by 


%u) * (63) 

which is now an explicit function 
of x. 


Sketch (i) indicates the varia- 
tion of the radius as a function of 
x/ Co for T 0 /mc 0 equal to 0.1 and 
values of B 0 /mc 0 equal to 0.5 and 
1 . 0 . 

Eelations (l4b) and ( 53 ) can be 
immediately employed for determining 
the pressure coefficient on the sur- 
face of the wing and on the body in 
region 2, because the portion of the body in that region is a circular 
cylinder. Use of these relations with equation ( 36 ), (57), and ( 5 8 ), then 
finally yields for the coefficient on the upper surface of the wing 

c p = m^rp^s - 3c 0 m) + 2(4s - c Q m + (c Q m -4s - 4y) + 



Sketch (i) 


c 0 m 


4s + 4 


R a 2 \ /s - c Q m -R g 2 /y 

7 )4 s - R 2 2 /y 


-4^' 


s - c Q m +R 2 2 /y 1 
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and for "the coefficient on the body 


£p 2 =m2T 


2(4s - 3Com) + 2(4s - c Q m) J + (c 0 m - 4s - 4fi 2 e i0 ) — ° m + R ^^ ~- + 

’ v a _ n _ c.iC7 


s - R 2 e 


/ 1 1 t_ -iSv -/ s - c n m + R 2 e -i ® , . , i0 N / a - c 0 m -R 2 e i0 

(c Q m - 4 b - 4R 2 e ) / 0 - ~ t TT 1 +(c Q m - 4s + 4R 2 e ) / + 

V s - R 2 e“ 1& v 8 R e i0 


( Co m - 4s + 4R 2 e- ie )y 


These equations are employed for 


> R ° . 
x ^ c o - m" + 


-i0 


(65) 


I Rp 2 ' CpOmr 

J 4 


the equality sign giving the value of x at the trailing-edge fuselage 
junction. The two formulas do not involve Mach number and apply for sub- 
sonic, transonic, and supersonic flight speeds. Comparison of formula 
(64) with equations (42b) and (45b) shows that for a fixed y in region 2 
the inequality 


°P > C Psub C Psup 


( 66 ) 


is satisfied when x > c 0 . Values calculated in region 2 for the subsonic 
and supersonic pressure coefficients along a section on the surface of the 
wing of the basic combination will thus be less than the values obtained 
for the coefficient Cp on the wing surface of the indented combination. 

In region 1, the body is not cylindrical but it is found that the 
squared terms in the pressure relation (37) may be neglected since they 
contribute only q uant ities involving the second and higher order in wing 

2 ^CD 

thickness that are small in comparison with - — • Even the first 

u o ox 

term in relation (37) gives arise, for the particular combination consid- 
ered, to some small terms of the second order. Such quantities, however, 
are also found to be negligible. The formulas to be presented here for 
the pressure coefficients in region 1 will therefore contain (like those 
for region 2) only terms of the same order as the thickness of the wing. 



Use of equations (l4b), ( 36 ), ( 53 ), (59), and ( 60 ) thus furnish the following final expres- 
sions for the pressure coefficients on the wing surface and on the body in region 1 of the combina- 
tion in sketch (h) : 

Spi = ^(ttCo - 4s + 4 Rl )/^I in - 8-7 (s - Rl )(c 0 m - s +R X ) + 


4(4s - 3c 0 m)tan“ ^ + F i(x,~y) - 2Fi(x,0) + Fj. (^x,~-)+F 1 (x,- ( 67 ) 


and 


~ w£t 

°Px " “ST 


-8 s/(s - Ri) (c 0 m - s +Ri) + 4(4s - Sc-mjtan ' 1 / — — — + F 1 (x,R 1 e' 1 ' 0 ) +F 1 (x,Rie' 

u V Co® - s + R X 


10 


) - 


SFiCx^O) +F 1 (x,-R 1 e 10 ) + Fj.(x, -Rxe -10 ) +Re 


2(= 0 m- ^ m 


e2i0 


U 219 -d 


where for c 0 m - s + 1 ] >0 

Fi(x,tj) - ( c 0 m - 4s + 4 tj ) J -° - m - 7 S -±l 

u V s - Tj 

inf 


10 


mC 0 (T) - Rl ) 


L2s ( c 0 m - s) + (2s - c 0 m) (tj + R l ) - 2R X + 2 -J (s - rj) (c Q ra- s + tj)(b - R 1 )(c 0 m - a + RJ J 
and for c Q m - s + tj < 0 


FjXxjTj) = (com - 4s + 4) 


J 


b - Com - TJ 


s - 


IT 2s(c 0 m- s) + (2 b - c 0 m)(T! +R X ) - 2 tjR x 

5" tan — f. ■ .. — 

2*7 (s - tj ) ( b - Com-Tj)(B -R 1 )(c 0 m- s + R x ) 


( 68 ) 


1 (69) 


<70) 


ro 

\o 
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NUMERICAL RESULTS AND DISCUSSION 


Graphs of the pressure coefficients for a wing alone and for the 
two wing -body combinations in sketches (f) and (h) are shown in figures 1, 
2, and 3« Figures 1 and 2 give plots for the coefficients along several 
spanwise stations of the basic combination (R 0 /c 0 m = 0.5) for subsonic 
and supersonic flow at m3 = 0.5, while figure 3 presents plots of the 
coefficient for the indented wing-body combination which are independent 
of Mach number. The plot in figure 4 gives the variation of the pressure 
on the body in the plane of the wing. 

Figures 1 and 2 show that, except at the wing-fuselage juncture, the 
pressure coefficients have a finite negative 6 value at the leading edge ' 
and increase to finite positive values at the trailing edge. On sections 
which are cut by the plane x = c Q (i.e., passing through the trailing- 
edge fuselage juncture) the slopes of the curves are discontinuous . It 
is also apparent that the effect of the presence of the body on the coef- 
ficient does not extend very far downstream beyond the trailing-edge 
fuselage juncture; the difference between the calculations for the wing 
on the combination and those for the wing alone, for instance, are too 
small to show up in the plots for sections more than one chord length 
from the body. 

•In figure 3j the pressure coefficients for nearly all sections on the 
wing of the indented combination are very close to the curve for the two- 
dimensional wing. The discontinuity in the slope of the curve that was 
quite noticeable along sections such as BB of figures 1 and 2 is far less 
apparent in figure 3 • 

The graph in figure 5 is the function which, according to equa- 
tions (27) , can be added to the values given in figures 1 and 2 for the 
subsonic or supersonic pressure coefficient along the section AA and BB 
of the wing to yield values of the coefficient at other Mach numbers. 
(Sections CC, DD, and EE lie in region 2 where the coefficients do not 

depend on Mach number.) For example, the pressure at the 

(T 0 / C o) m 

leading edge along BB in subsonic flow for mp = 0.4 is. using figures 1 
and 5 in conjunction with equations (27), -4.1 -0.2 = -4.3- 

Figure 6 shows isobaric charts of the pressure coefficients for the 
wing without body in subso nic and supersonic flow at mp = 0.5^ also 

6 The occurrence of a negative pressure at the nose of a Joukowski- 
like section is the result of the thin-airfoil simplification. A more 
accurate theory would show a small region of positive pressure (i.e., 
a st agnat ion point with the maximum value equal to the impact pressure 
of the component of stream velocity normal to the edge) . 
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a chart giving the -two-dimensional results. The figures 6(a) and 6(b) 
indicate that the pattern of the isobars in the (region behind the trailing- 
edge fuselage juncture is essentially the same as the tiro-dimensional in 
figure 6(c), but that in the region near and upstream of the juncture, a 
marked deviation from the straight isobars of the tiro -dimensional case is 
evident . 

Isobaric maps are also shown in figure 7 for the basic wing-body 
combination in subsonic and supersonic flow, and for- the indented wing- 
body combination. Figures 7(a) and 7(b) illustrate that in the region • 
downstream from the trailing-edge fuselage juncture (even near the juncture 
itself) the isobars are not much different from those for the wing alone, 
and that the body therefore has little effect in this region. In the 
region adjacent to the body, the pattern of the isobars is qualitatively 
similar to the case for the wing alone but the pressures are lower. A 
remarkable difference between the chart in figure 7(°) for the indented 
combination and those in figures 7(a) and 7 (b) for .the basic combination 
is that the isobaric pattern on the wing for the indented combination is 
essentially two-dimensional over practically the whole wing. 

Examination of the three cha r ts in figure 7 also indicate that the 
maximum negative pressure on the wing occurs at the leading edge near the 
boundary between regions 1 and 2. In fact, it can be shown that the max- 
imum occurs at the boundary. In view of the inequality (66), indentation 
in accordance with the area rule reduces the maximum perturbation veloci- 
ties on the wing. 


Ames Aeronautical Laboratory 

National Advisory Committee for Aeronautics 
Moffett Field, Calif., Oct. 7, 1954 
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Figure 1.- Subsonic pressure coefficient on upper surface of wing of the 

basic combination; m3 = 0 . 5 . 






















































